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Spin polarization effects in micro black hole evaporation 
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We consider the evaporation of rotating micro black holes produced in highly energetic particle 
collisions, taking into account the polarization due to the coupling between the spin of the emitted 
particles and the angular momentum of the black hole. The effect of rotation shows up in the helicity 
dependent angular distribution significantly. By using this effect, there is a possibility to determine 
the axis of rotation for each black hole formed, suggesting a way to improve the statistics. Deviation 
from thermal spectrum is also a signature of rotation. This deviation is due to the fact that rapidly 
rotating holes have an effective temperature T e g significantly higher than the Hawking temperature 
Th- The deformation of the spectral shape becomes evident only for very rapidly rotating cases. 
We show that, since the spectrum follows a blackbody profile with an effective temperature, it is 
difficult to determine both the number of extra-dimensions and the rotation parameter from the 
energy spectrum alone. We argue that the helicity dependent angular distribution may provide a 
way to resolve this degeneracy. We illustrate the above results for the case of fermions. 



Introduction. Within the context of TeV-scale grav- 
ity [H, 0, H] , the possibility that colliders or cosmic ray 
facilities may observe micro black holes has attracted 
enormous attention [f| 0, 0, H, H|. A close look at the 
limits on the fundamental Planck scale shows that a win- 
dow of about 5 TeV is still open for the LHC to observe 
such exotic events floj |. while the window is much wider 
for cosmic rays. Micro black holes with even higher en- 
ergies could be produced from the collision of a cosmic 
ray with an atmospheric nucleon, a dark matter particle, 
or another cosmic ray (Ref. HH gives some up-to-date 
estimates). 

In this paper, we consider micro black holes result- 
ing from the collision of two particles at energies much 
higher than the higher dimensional Planck mass Mp. We 
have in mind models with Mp of order of a few TeV 
and the standard model confined on a 3-brane, embed- 
ded in a (4 + n)-dimensional bulk. These black holes 
have horizon radius smaller than the size of the extra di- 
mensions, and arc expected to follow balding, spin-down, 
Schwarzschild, and Planck phases. Micro black hole for- 
mation has been studied both analytically [l2j and nu- 
merically [H, [l4j, and their evaporation has also been 
the subject of considerable attention (See Refs. [l(| EH 
for review). Previous work suggests that micro black 
holes will emit mostly brane modes [l?], tHI, and the 
deviations from the blackbody spectrum have been in- 
vestigated using numerical and semi-analytical meth- 
ods pi m mi, m, m nnm. 

Radiation. In this paper we analyze the fermion emis- 
sion from spinning evaporating micro black holes. As- 
suming that the black hole horizon is significantly smaller 
than the extra dimensions, we approximate it by a vac- 
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uum higher dimensional Kerr [2j 
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where A = r 2 +a 2 — Mr 1 ~ n and E = r 2 + a 2 cos 2 8. Mp is 
normalized to one. Since we are interested in the visible 
brane modes, the background spacetime will be given by 
the projection of the above metric on the brane. 

Massless fermions emitted by the black hole are de- 
scribed by the Dirac equation: 

where ip is the Dirac spinor wave function, a set of 
tetrads, the spin-affinc connections determined by 
T M = 7 a 7 b ^a fcAI /4 , with Uabti being the Ricci rotation 



are chosen to sat- 
g^ v , with being the 
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coefficients. The matrices 7^ 
isfy the relation 7^7^ + 7 I/ 7 A ' 
metric on the brane. 

The Dirac equation for massless fermions on a Kerr 
background has been studied extensively in four and 
higher dimensions @, [H, HE IH HI SO, EH]. Here we 
briefly describe the calculation, closely following the ap- 
proach of Ref. 28]. Due to the symmetries of the Kerr 
spacetime, the spinor wave function factorizes as 
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where the + and — signs refer to negative and positive 
hclicitics, respectively. We illustrate the results for the 
case of negative helicity. The positive helicity case can be 
obtained by a trivial chirality transformation. The field 
d> takes the form 
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FIG. 1: The normalized energy spectrum of the emitted 
fermions. The horizontal axis is rescaled by the effective tem- 
perature determined by fitting the data by a black body pro- 
file. The overall amplitude is also normalized since the abso- 
lute magnitude is not observable. The upper and the lower 
panels are the plots for n = 2 and n = 4, respectively. The 
rotation parameter a is set to 30%, 50% and 70% of a cr u 
(left) and a max (right). 

and the normalization factor is 9i~ x = A 1 / 4 (r + 
ia cos 0) 1 / 2 sin 1 ' 2 9. The angular and radial modes obey 

(±± U asine T ^)s T {e) = ±kS ± (6), 

T A~ + fl2) " m&) ) i?T(r) = ^A^R^r) , 

where n is a separation constant. Supplemented with 
regularity conditions at 6 = and 7r, the set of angular 
equations provides an eigenvalue problem, which deter- 
mines k [30(. In order to compute the particle flux, we 
need the solution to the radial equation supplemented by 
ingoing boundary conditions at the horizon 

J?_ ~ , R+ ~ e-™ T * , for r* -> -oo 

where a) = cj — maj(r\ + a 2 ), and r* is defined by 
dr*/dr = (r 2 + a 2 )/ A. The number of particles emit- 
ted, for fixed frequency co, is distributed according to the 
Hawking radiation formula. For negative hclicity modes, 
the angular distribution reads: 

dN = 1 V^io (n\\2 gVn n x 
dcudcose 27rsin6»^ 1 ~ l ;l e "/ T H + 1 ' ^ ' 

where TW = _!_ ("+ 1 ) r h+(»- 1 ) a j s fo e Hawking tem- 

perature, and the grey-body factor er ; is the squared 
amplitude of the transmission coefficient of an incoming 
wave from r = oo (see Ref. (3l|). 



Results. The initial angular momentum of the pro- 
duced black holes J = 2aM/(n + 2) is restricted by re- 
quiring the impact parameter b = J/M to be smaller 
than the horizon radius r^, determined by A(r/j) = 0. 
Then, the maximum value of the rotation parameter a 
turns out to be a max = ^y^r/j 20]. The upper bound 
on J might be even lower for n > 2. In fact, there ex- 
ists a critical value for a, a cr u = (n + l)(n — l) _1 r 2 , 
where \d(T,Qh)/d(M,J)\ vanishes. If the same argu- 
ment as in the case of black branes applies, black holes 
with a > a cr it suffer from the Gregory Laflamme insta- 
bility (See also Ref. [32l]). Then, a cr a represents the max- 
imal value below which the higher dimensional Kerr so- 
lution is adequate. The value at which the dynamical 
instability is expected to set in may be slightly differ- 
ent from a cr it, which only represents an indicative esti- 
mate of the critical value. Interestingly a cr n < a max (for 
n = 2, 3, 4 extra dimensions, a cr a = 1.09, 1.07, 1.06, 
whereas a max = 1.25, 1.89, 2.46). Although it is widely 
believed that a dynamical instability exists, the value of 
a crit obtained above is only heuristic. Thus, we con- 
sider two possible cases: the maximal value allowed for 
a is a cr it or a max . A set of representative values for the 
parameter a is chosen as a/a max — 0.3, 0.5, 0.7, and 
a/a crit = 0.3, 0.5, 0.7. M is set to unity. Another 
natural choice to present the results would be to use the 
impact parameter b. In our case, we use the ratios a/a cr it 
and a/a max , which correspond respectively to b/b cr it and 
b/b max , where b cr n and b max are the values for the im- 
pact paramter corresponding to the critical and maximal 
cases respectively. 

Having fixed a in the above way, we compute the en- 
ergy spectrum, shown in Fig. [1] We normalize the hori- 
zontal axis by using an effective temperature T c g deter- 
mined by fitting the data by a blackbody spectrum pro- 
file. The effective temperature T e g is much higher than 
the Hawking temperature as shown in Fig. However, 
the spectral shape is not so different from the thermal 
one except for the cases with a « a max (Fig.QJ. In previ- 
ous work, [25|, [2(| , the enhancement of emission at large 
frequencies is reported. However, the deviation from a 
blackbody spectrum was not quantified. 

We find that the renormalized spectra are enhanced 
for both lower and higher frequencies compared with the 
black body spectrum at T = T e g (thick line). Except for 
very large values of a, we have shown that the obtained 
spectra can be fit well by superpositions of black body 
profiles with width of about 2tt H x T e g. The intuitive 
reason for the enhancement of the effective temperature 
is that the motion of the hypothetical emitting surface on 
the rotating black hole, relative to observers at infinity, 
causes an additional blucshift factor (which varies from 
place to place) . This is because corotating emitted parti- 
cles encounter less suppression from the statistical factor. 
This can be made precise by closely inspecting the com- 
bined behavior of the greybody and statistical factors. 
The dominant contribution to the spectrum comes from 
the I = m modes and for larger values of the rotation 
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FIG. 2: The effective temperature normalized by the temper- 
ature at a — versus the rotation parameter a normalized 
by a cr it (left) and amax (right). The curves represent the 
Hawking temperature. 



the contribution to the spectrum from such modes, with 
large I, is non negligible. 

However, because of the change in the temperature 
and the rotation parameter during the evaporation, the 
broadening of the spectrum due to the rotation will not 
be identified straightforwardly. Wiggles can be also seen 
in the spectrum for a small number of extra-dimensions. 
However, wiggles are likely to disappear as T and a 
change during the evaporation. 

When the rotation velocity is high, the deviation from 
the thermal spectrum is much clearer. As a novel signa- 
ture, we find that the spectrum is sharply cut off at high- 
frequencies for rapid rotation. This new signature may 
survive even after we take into account the superposition 
of spectra along the evolutionary track of an evaporating 
micro black hole. This highly spinning regime is realized 
for a > a cr i t . 

In Fig. [31 the angular distribution of negative helicity 
particles is displayed for various parameters, setting uj 
to a representative frequency uj. The value uj is chosen 
by requiring that the fraction of particles emitted with 
frequency below uj, N(uj) = J Q dN, to be 0.5. 

The emission turns out to be suppressed in the direc- 
tion anti-parallel to the black hole angular momentum. 
For rapid rotation, the particles tend to be emitted to- 
wards the equatorial plane. This concentration in the 
rapidly rotating case can also be seen in the helicity inde- 
pendent angular distribution [26j |. The emission around 
both poles looks suppressed, but the observed apparent 
suppression is simply caused by the large enhancement of 
emission in the directions close to the equatorial plane. 
The asymmetry in the helicity dependent angular distri- 
bution is visible even for relatively slow rotation and be- 
comes evident as a increases. Note that, for very fast ro- 
tation, the concentration of the emitted particles toward 
the equatorial plane, observed in the angular distribu- 
tion, may affect the features of cosmic ray air showers 
mediated by black holes. 

For slow (rapid) rotation, the asymmetry decreases (in- 
creases) as the number of extra-dimensions n grows. This 
tendency may be used as an indicator to discriminate 
scenarios with different number of extra-dimensions. It 
is also worth mentioning that for a/a cr it fixed the peak 



n/a„, 0.3 0.5 0.7 a/acnt 0.3 0.5 0.7 

n = 2 18.20 15.17 9.47 n = 2 20.68 16.20 13.17 

n = 3 19.93 13.43 8.19 n = 3 25.47 19.32 15.34 

n = 4 20.03 10.97 7.50 n = 4 29.84 21.75 17.14 



TABLE I: Estimate of 5 in degrees for the curves of fig. [3] 

position of the helicity dependent angular distribution is 
almost independent of n as shown in upper panel, Fig. [3] 

Statistics. If we can align the direction of the axis of 
rotation of the black hole for various events even approx- 
imately, we can collectively use the experimental data 
to achieve high statistics for the angular distribution of 
emitted particles. The LHC may allow to perform such 
measurements. 

In the following, we will provide an estimation of the 
error in the determination of the axis of rotation, as- 
suming that N particles are emitted per black hole. Let 
P(f2) be the angular distribution of the emitted parti- 
cles. We expand it in terms of Legendre polynomials as 
P(f2) = ^2 CiPi (cos 6), and consider the identification 
of the direction of angular momentum based on I = 1 
(dipolc) and I = 2 (quadrupolc) moments. For the dipole 
and quadrupolc estimators, respectively, the errors in 
the estimated direction, Sd and S q , can be evaluated as 

S d = ^cf and 5 l = jf§^$> with C = Co/3 + 2C 2 /15. 
Combining the dipolc and quadrupole estimators, the la- 

error in total can be reduced to 5 = (S^ 2 + 8~ 2 ) 
Assuming that the angular distribution shown in Fig. O 
the error S in degrees for N = 100 is summarized in Ta- 
ble UJ Here we restricted our consideration to the dipole 
and quadrupolc estimators, but more sophisticated sta- 
tistical analyses may reduce the error. 

Discussions. In the collision of two particles at trans- 
plankian energy, a rotating black hole is expected to form 
and decay. We studied possible signatures of rotation of 
such black holes in departure of the energy spectra from 
the thermal profile, and in the features of the helicity 
dependent angular distributions. 

Black hole thermodynamics seems to suggest the pres- 
ence of an instability for a > a cr i t . This critical value 
a cr it is smaller than the maximal value a max allowed by 
the kinematical requirement of formation of a black hole 
in the collision of two particles. As far as a < a cr a, the 
shape of the energy spectrum is almost independent of 
n. The largest dependence on n will appear in the ef- 
fective temperature as shown in Figs. [21 However, this 
n-dependence must be interpreted with caution. When 
the ratio a/a cr it is fixed, the enhancement of the effec- 
tive temperature is larger for a smaller number of extra- 
dimensions. While the tendency is completely opposite 
if the ratio aja max is fixed. Hence, under the situation 
in which the true maximum value of a is unknown, it is 
rather difficult to extract the information about the num- 
ber of extra-dimensions without changing the energy of 
colliding particles @ . 
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FIG. 3: The angular distribution of emitted negative helicity fermions. In each plot the cases with n — 2,3 and 4 are shown 
simultaneously. The reference values of the rotation parameter are a cr it and a max in the upper and lower panels, respectively. 
The rotation parameter a is set (left to right) to 30%, 50%, 70% of the reference value. 



The peak position of the helicity dependent angular 
distribution may give a valuable information, because it 
seems to be a good indicator of a/a cr i t (or a itself since 
a cr it is always close to 1). Moreover, the amplitude of the 
anisotropy depends on the number of extra-dimensions. 
Hence, measuring the helicity dependent angular distri- 
bution may provide a very important signature to extract 
the value of n. To develop analysis of this kind based on 
experiment, we need to coherently accumulate data from 
many events. For this purpose, it is necessary to iden- 
tify the rotation axis of the formed black hole for each 
event. We have demonstrated that this identification is 



marginally possible if we can detect 0(100) particles. 
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